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1. Introduction 

Spintronics in molecular conductors is a field attracting more and more attention, 
both from fundamental physics as well as from application-oriented material science 
PQ. Here the quantum- mechanical electronic spin is the central object controlling 
transport properties. For a conductor sandwiched between ferromagnetic leads, a 
different resistance can be observed depending on the relative orientation of the lead 
magnetizations. Quite often, the resistance is larger in the antiparallel configuration 
than in the parallel one, but sometimes also the reverse situation can be observed. It 
is useful to define the tunnel magnetoresistance (TMR), p t = (Rap — Rp)/Rp, as the 
relative difference between the corresponding resistances. 

A particularly interesting material in that context is provided by carbon nanotubes 
(CNTs), see Refs. [21 Ej for general reviews. Quite a number of experimental studies 
concerning spin transport through individual multi- (MWNT) or single-walled (SWNT) 
nanotubes contacted by ferromagnetic leads have been reported over the past few years 
jH El El HI El EE UH] • In particular, the experiments of the Basel group [TDJ use thin- 
film PdNi alloys as ferromagnetic leads in order to contact either SWNTs or MWNTs, 
where the shape anisotropy and the geometry of the setup allow for the study of the spin- 
dependence of electrical transport. These experiments have revealed oscillatory behavior 
of the TMR as a function of the external gate voltage. Similar oscillations were predicted 
as a consequence of the gate- voltage-tunable Rashba spin-orbit (SO) interaction [TT| IT2] 
in a classic paper by Datta and Das some time ago [13] . Since Datta-Das oscillations 
have still not been observed experimentally so far, a thorough theoretical investigation 
of this effect in nanotubes is called for and provided here. Unfortunately, from our 
analysis below, we find that the weakness of SO couplings in nanotubes excludes an 
interpretation of these data in terms of the Datta-Das effect - they can, however, be 
explained in terms of quantum interference effects ^0]- Nevertheless, we show that the 
presence of multiple bands in CNTs is not detrimental, and under certain circumstances, 
the effect may be sufficiently enhanced to be observable, e.g., by a tuning of the number 
of bands via external gates along the lines of Ref . . In the original Datta-Das proposal 
[T3] . subband mixing was ignored so that different channels just add up coherently, but 
subband mixing has later been argued to spoil the effect [TH1 [THJ - In CNTs, the special 
band structure requires a careful re-examination of the Datta-Das idea in this context, 
and we shall show that the arguments of Refs. [T5] ITfi] do not necessarily apply here. 

Recent theoretical studies of spin-dependent transport in CNTs have mainly focused 
on the single-channel limit, taking into account electron-electron interactions within the 
framework of the Luttinger liquid theory [TH [TBI QUI EOl E5 (see also [221 E31 E3 for 
related discussions on interacting quantum wires with Rashba SO coupling). Here we 
confine ourselves to the noninteracting problem in order to not overly complicate the 
analysis, but study the many-band case and details of the band structure. Interactions 
can be taken into account within the Luttinger liquid approach at a later stage, and 
may enhance the effect of SO couplings [221 123- We shall also neglect disorder effects. 
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Mean free paths in high-quality SWNTs typically exceed 1/zm, while in MWNTs this 
may be a more severe approximation for some samples. However, high-quality MWNTs 
with ultra-long mean free paths have also been reported recently 

The structure of this paper is as follows. In Sec. |2 we derive the Rashba spin- 
orbit hamiltonian from microscopic considerations. The resulting tight-binding SO 
hamiltonian will be studied at low energy scales in Sec. El where we derive its continuum 
form. In Sec. 01 the consequences with regard to Datta-Das oscillations in the TMR 
are analyzed. We shall always consider the zero-temperature limit, and (in most of the 
paper) put h = 1. 

2. Rashba spin-orbit coupling in nanotubes 

We start by noting that transport effectively proceeds through the outermost shell of 
a MWNT only, such that we can take a single-shell model even when dealing with a 
MWNT. Experimentally and theoretically, it is understood that such a model works 
very well in good-quality MWNTs [2], essentially because only the outermost shell 
is electrically contacted and tunneling between different shells is largely suppressed 
[2Z1I2H|- Naturally, a single-shell description is also appropriate for SWNTs, where we 
assume a sufficiently large radius R such that occupation of multiple subbands can be 
possible. (For a MWNT, R denotes the radius of the outermost shell.) Depending 
on the electrochemical potential \x (doping level), we then have to deal with N spin- 
degenerate bands. We assume full quantum coherence (no dephasing), so that the usual 
Landauer-Btittiker approach applies, and exclude external magnetic fields or electric 
field inhomogeneities, say, due to the electrodes. We proceed to derive the Rashba SO 
interaction, H so , for this problem. Notice that this is different from the intrinsic atomic 
SO interaction discussed in Refs. (TBI EH]- In particular, the SO coupling in Refs. 
vanishes in the limit of large radius, which is not the case for the Rashba SO coupling we 
discuss below. Though Ando's SO coupling [TH] could straightforwardly be included in 
our analysis, being gate- voltage independent it could not change our conclusions relative 
to the gate-voltage dependent oscillations in the magnetoresistance and is neglected in 
what follows. 

We first define a fixed reference frame S = {Y, Z, X}, with unit vector X pointing 
in the axis direction and Z perpendicular to the substrate on which the CNT is supposed 
to be located. Next we introduce a second, local reference frame «Sj = {pi, £j, X} relative 
to each lattice site Ri on the tube surface, where pi and i% are unit vectors along the 
local normal and tangential (around the circumference) directions at Ri, respectively. 
Using polar coordinates in the plane transverse to the tube axis, the relation between 
S and Si is given by 

Pi = cos ipiY + simpiZ, £j — — sin cpiY + cos (piZ. (1) 

The position vector of a given carbon atom can then be written as Ri = Rpi + XjJf . 
For later convenience, we introduce also another reference frame. For each pair of sites 
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Ri and Rj, we define 

Rij = Ri — Rj = X^X + pij, (2) 

and denote the direction perpendicular to and X as pjj. Then {pjj, Pij, X} constitutes 
a new local frame «Sy, and one has 

p^ = cos[(^ + <fj)/2] Y + sin[(^ + <pj)/2] Z, (3) 

Pij = - sin[(v?; + ipj)/2] Y + cos[0; + <fj)/2] Z. 

The 2p z orbital at position Ri can then be represented as 

Xi (f- Ri) = a(f- Ri) ■ pie-W-^, (4) 

where 4a = (27rao)~ 1//2 , f3 = (2a ) _1 , a = h 2 /me 2 = 0.53A is the Bohr radius, and m is 
the electron's mass. We introduce an index i on the orbital in order to keep track of the 
atom at which it is centered. The wavefunction (jlj is expected to be highly accurate 
for not too small R, where hybridization with the sp 2 orbitals is negligible. 

At large distances from the tube, external gates generally produce an electric field 
perpendicular to the tube axis and the substrate. As it has been shown in detail in 
previous works [3U1 Ell > polarization effects of the CNT itself due to a transverse field 
result in a reduction of the externally applied field described by 

1 

l + 2c7 0yy /R 2 ' 

where a 0yy is the unscreened transverse static polarizability. Since a 0yy is approximately 
proportional to R 2 , the factor in front of E ext practically equals a constant, ~ 0.2 [30J. 
Then, assuming homogeneity, the electric field due to the gate can be written as 

E = E Z, (5) 

which in turn produces the (first-quantized) Rashba spin-orbit interaction JH^j. With 
standard Pauli matrices a acting in spin space, 

eTi -* 

Hso = ^ E - {&X ^ (6) 
We proceed to derive the second-quantized spin-orbit hamiltonian within the tight- 
binding approximation. For that purpose, we need the matrix element of the momentum 
operator between two 2p z orbitals p^ = (Xi\p\Xj)i from which we get the following form 
for the SO hamiltonian: 

H so = gJ2 c l [(<? x Pij) ■ Z] Cj, (7) 

ij 

where the fermionic operator Ci a destroys an electron with spin a =f, J. in the 2p z orbital 
centered at Ri, and g = Eo/4m 2 c 2 . For calculational convenience, the matrix element 
Pij can be written as gpij = i(vij + Utj), where the spin-orbit vectors Vij and Uij are 
defined as 



Eq — - — jr^E ext , 



v 



u 



-gaj d s f Xl (r- R^p^-^, (8) 

g(3 [ d 3 r Xt (f- R^l^-yjif- Rj), (9) 
J \r-Rj\ 
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Note that the modulus of and Uy has dimension of energy, and their sum (but not 
necessarily each term separately) is antisymmetric under exchange of % and j. 

We first observe that the spin-orbit vectors connecting a site with itself clearly 
vanish, since (Xi\p\Xi) = 0. Let us then discuss spin-orbit vectors connecting different 
sites. Since the orbitals (jlj) decay exponentially, it is sufficient to consider only the case 
of nearest neighbors. We start with Vy. Shifting r — > s + Ri in Eq. (jSJ) and using Eq. (jlj), 
we obtain 

2 n. f ^ 3 ?r?. a,\ P -P' e -P\***i\ 



^ = -ga pj j d s(s- pi)e 
Using s — s\\Rij + s±, we then rewrite the above integral as 
J d 3 s( Sll R tJ + s ± ) ■ p,e-^e-V( s ii+ d ) 2 +4 ) 

where we use \Rij\ = d, with the nearest-neighbor distance among carbon atoms in 
graphene d = 1.42 A. Note that (3d = 1.34. The second term in the brackets is odd in 
s± and thus vanishes, and we obtain 

v i: j = -ga pj— sin ( )d 7 , (10) 



(/ v 2 
f sin 2 (^ 



where we have used R^ ■ Pi = sin 2 ( y ' 2 yj ) and the dimensionless numerical factor 7 



7o = j dxd ! /^K^v / W e ^v / MW. 
For Vji, we find 

Vji = ga pi— sin ( — -)d 70 . 



Notice that, up to higher orders in d/R, the unit vectors p^j can be replaced by p^-, 
which makes clear that v^ is normal to the tube surface. Now | sin[(yjj — (pj)/2]\ varies 
between zero (when the two sites are aligned in the axis direction) and d/2R 1 (when 
the two sites are aligned in the circumferential direction). Thus, to zeroth order in d/R, 
v^ vanishes: it is a pure curvature effect, peculiar of nanotubes, which does not exist in 
graphene. In practice, v^ is tiny and certainly subleading to Uij, which turns out to be 
of order (d/R)°. We shall therefore neglect it in what follows. 

Let us now turn to Uij. We shift r — ► s + (Ri + Rj)/2 in Eq. 0, and rewrite as 
the sum of two terms: 

uf =g(3 f d 3 s Xi (s- R %3 j2) Xj (s + R^/2) | (11) 
J \s + Rij/2\ 

uf = ^4 / d?sUs- 4/2) Xj (s + Rij/2) 1 (12) 
2 J \s + Rij/2\ 

Writing again s = s\\Rij + s±, the computation of the above integrals leads, to the lowest 
non- vanishing order in d/R, to the following expressions: 



u\j> = gf3a 2 R ij d 4 'y 1 = Wi-R^-, (13) 



t\j = gj3a 2 Rijd ir 

Uif = ^® 2 Rijd A l2 = u 2 Rij, (14) 
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with the dimensionless numerical factors 



xz 2 e -Pdy/{x-l/2)l+yZ+zZ e -fldy/{x+l/2)2+yl+zl 

7i = / dxdydz , = (15) 

y/( X + 1/2)2 + y 2 + Z 2 

0.0375, 



and 



z 2 e -f3d^/{x-l/2)lW+z 2 e -Pd^/{x+l/2)lW+z 2 

72 = / dxdydz = (16) 

J ^(x + 1/2) 2 + y 2 + z 2 

~ 0.3748. 

To lowest order in d/R, it does not make a difference whether we take the tangent unit 
vector at Rj, Ri, or at (Ri + Rj)/2. Hence we may write pij — > e v , where e v is the unit 
tangent vector at (Ri + Rj)/2. We then get SO couplings along the axial and along the 
circumferential direction, 



Uij = u 



(Ri, ■ X)X + (Rij ■ e v )e v J , (17) 

with u = U\ + «2- Note that we have neglected a tiny component of R^ normal to 
the tube surface. The above discussion then results in the tight-binding hamiltonian 
H = H + H so , where 

r,a 

with t fa 2.7 eV Here the f denote all sublattice-A tight-binding sites of the lattice. 
Furthermore, the 5 a =i,2,3 are vectors connecting r with the three nearest-neighbor sites 
which are all located on sublattice B jH]. Since we consider the limit d/R 1, the S a at 
each site effectively lie in the tangent plane to the tube surface at that site. The Rashba 
spin-orbit hamiltonian then reads 

H so = iu J B ~ +Sa [(ff x [(5 a ■ X)X + (5 a ■ e v )e v ]) ■ z] c A> ? + h.c. (18) 

r,a 

3. Continuum limit 

Since we are interested in the low-energy long-wavelength properties, we now expand 
the electron operator around the Fermi points K, K' in terms of Bloch waves 



J= = e iK - ? F lp (r) + e- iK - r F 2p (r), (19) 

where S = \/3a 2 /2 is the area of the unit cell, a = \/3d, and p = A/B is the sublattice 
index. The F ap are slowly varying electron field operators, and we choose the Fermi 
points at K = (An/3a, 0) and K' = —K 0. We then expand F(f+6) ~ F(f) + 5-VF(r) 
and use the bond vectors 

5 = ^(0,-1), & = |(1,1/V3), ^8 = |(-1,1/V3). (20) 

These vectors are given in a fixed reference frame for a 2D graphene sheet, and we then 
must perform a rotation to longitudinal and circumferential directions via the chiral 
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angle. This rotation results in fixed phases that can be absorbed in the definition of F ap 
and do not appear in final results. This is of course expected from the U(l) symmetry 
emerging at low energies in the dispersion relation of graphene |3j. After some algebra, 
the usual Dirac hamiltonian for the kinetic term follows, 

H = vJ d 2 fF^ [(T ® r 2 ® a )(-id x ) + (T 3 ® n <g> a Q )(-id y )\ F, (21) 

where v = y/3at/2 ~ 8 x 10 5 m/sec is the Fermi velocity, and x,y are longitudinal and 
circumferential coordinates, respectively, with < y < 2irR. Finally, Tj and Xj are also 
Pauli matrices that now act in the space of Fermi (K, K') points and sublattice space 
(A, B), respectively. For i — 0, these are defined as 2 x 2 unit matrices. 

The low-energy limit of the SO term (|18j) can be obtained in the following way. 
First we observe that 

'(<?„ • X)X + (5 a ■ e v )e v ] x Z = -(5 a ■ X)Y - sm(y/R)(5 a ■ e v )X. 

Here the only approximation is the assumption that the bond vectors lie in the plane 
tangent to the nanotube surface at r. Second, by using the bond vectors (j2~Uj) and taking 
into account the chiral angle 77 between the fixed direction on the graphite sheet and 
the circumferential direction on the nanotube, one obtains 

£ C ks a • *) c ? « 4^ (Fke- iT >F M + Fl 2 e^F A2 ) , 

a 
a 

Notice that we take into account exactly the relative orientation of the bond vectors 
with respect to the directions X and e v for a generic nanotube, which is encoded in the 
chiral angle 77. The constant phases e ±lv can be absorbed by appropriately redefining 
the operators as F A2 — > e~ lv F A 2 and F S1 — > e~ ir, F B i, and the final result can be written 
down in the form 



U||T ® n ® (7 2 + ui sin(y/i2)T 3 ® r 2 ® <ri T, (22) 



with My = mj_ = 3du/2. For the sake of generality, we continue to use different coupling 
constants u± and u\\. It is worthwhile to mention that the leading term for the Rashba 
spin-orbit coupling in a CNT, Eq. (|22|). does not depend on longitudinal momentum. 
This is due to the peculiar band structure of graphene with its isolated Fermi (K) 
points. In the above derivation, we also find terms that are linear in momentum, 
i.e., contain spatial derivatives of the electron operators. Such terms only produce 
tiny renormalizations of the velocities and will be neglected here. The second term in 
Eq. allows for spin flips and mixes transverse subbands. 

From now on, for simplicity, we consider just a single Fermi point, say, K. After 
the global SU{2) rotation o\ — > cr 2 — » cr 3 in spin space, we get in compact notation 

Ho =v [-indy - ir 2 d x ] . (23) 
*Hso = M||Ticr 3 + u± sm(y/R)T 2 a 2 . (24) 
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Note that the exact spectrum of Ti.Q+Ti. so with u± = can be obtained straightforwardly. 
In general, however, due to the smallness of the SO coupling (see below), it is enough to 
treat TL so perturbatively. The following detailed derivation is then necessary to correctly 
evaluate the effect of the SO coupling, and moreover it is interesting and important for 
the generalization to the interacting case, and for the analysis of features involving the 
electron wavefunction (as for instance electron-phonon interactions). 
The eigenvalues of 7io are given by 

e a na(q) = av^Jk 2 ± (n) + q 2 = ae n (q), (25) 

where k±(n) = (n+n )/R denotes the transverse momentum, q the longitudinal one, a = 
± labels the conduction/valence band, and a = ± the spin. Here uq = for intrinsically 
metallic shells, but generally it can be taken as < no < 1/2 to take into account 
chirality gaps or orbital magnetic fields along X. The transverse subbands are labeled by 
integer values n = 0, 1, 2, . . . , W-l, where Af = 2(N 2 +M 2 +NM)/gcd(2M+N, 2N+M) 
for (N, M) tubes j^j. M is typically much larger than the actual number iV = [k F R] of 
occupied subbands, where we define kp = fi/v with the doping level fi that we assume 
positive here. The velocity v n for electrons in subband n at the Fermi level (in the 
absence of H so ) and the corresponding Fermi momentum q n are then given by 

v n = v\]l - [{n + n )/(k F R)] 2 , q n = k F v n /v. (26) 

The eigenvalues (J23j) are spin- independent and thus doubly degenerate. The 
corresponding eigenstates are denoted \nqao~), where \n) and |g) are respectively plane 
waves in circumferential and longitudinal direction. In coordinate representation they 
read 

e ik ± (n)y 

*Pn Q aa(x, y) = (x, y\nqaa} = -j==e iqx £ na (q) <g> x*, (27) 
with the bispinor (in sublattice space) 

e eW 1 e ^)- lMz^ (28) 

U,a=±[q) - ^ y ±e -ie n{q) ,2 ), e - €n{q) . (2»j 

A different, and here more convenient basis is given by the sublattice states \nqpa). 
Their coordinate representation is 

ik ± (n)y 

e iqx t P ®X«i (29) 




V2^R 



B 




Their usefulness stems from the fact that the \nqpo~) can be factorized as 

\nqpa) = \n)\q) ® \pa), \pa) = £ p ® Xa, (30) 
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where \pa) is independent of n and q. Using this basis, we can expand the field operator 
F(r) on the tube surface as 

F (r) = J2 [ p-^n qP a{x,y)c npu {q) =^2F n (x)(y\n), (31) 

where the operator c npa (q) destroys an electron in the state \nqpa), and we introduce 
the ID field operators F n {x). Alternatively, using the basis of eigenstates of H , F(r) 
can be expanded as 

F ( f ) = [ P-' l PnqacT(x,y)c naa (q), (32) 

where the operators c naa (q) destroy conduction (a = +) or valence (a = — ) electrons 
with spin a in subband n. Notice that in what follows the spin index is left implicit. 
The relation between the operators c na and c np is easily found to be 

/ c n+ {q) \ J_ / e-^/ 2 e^/ 2 \ ( c nA (q) 
V c n -(q) ) V2 { e-^ 2 -e^ 2 ) \ c nB {q) 

We now proceed by treating the spin-orbit hamiltonian using perturbation theory. 
First, we diagonalize H — /j,N for a fixed transverse subband n, 

- /iiV» =v J dxFl[k ± (n)n + {-id x )r 2 - p]F n 
dq 



(33) 



J2 / ^K>(<?) -tA C L C r, 



a=± 

Next we expand around the Fermi points ±g n defined in Eq. (|2fij) . which introduces 
right- and left-movers, r = ± = R/L, as the relevant low-energy degrees of freedom. 
For small deviations k from ±q n , Taylor expansion yields e n (±q n + k) ~ /i ± v n k, where 
v n is given in Eq. ()26|) . Since we assumed /i > 0, we may now restrict ourselves to the 
conduction band, a = +. For the hamiltonian, we then obtain 

dk 



- /iiV^ = £ v n f -(rk)cUk)c nr (k) 
J2v n Jdx ipl r (-ird x )if; nr , 



r=± 



where c nr (k) = c n+ {rq n + k) and 4> nr (x) = J #e c nr {k). This introduces i?/L-moving 
ID fermion operators for each subband n (and spin a). The relation of these ID fermions 
with the original operator F n (x) is given by 

^ e ikx / e i0 n (q n )/2 

2^71 V e" 

dke ikx ( e-^" 1 """- 
2^72 V e * 

Notice that, while in general the unitary transformation from sublattice space to 
the conduction/valence band description depends on longitudinal momentum, in the 
continuum limit, one can use the transformation directly at the Fermi momenta. This 



F n (x) = j -— [ _ WM/2 c nR (k) 



+ Mqn)/2 c nL (k). (34) 
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is consistent with the neglect of band curvature effects implicit in the linearization of the 
dispersion relation, which is unproblematic away from van Hove singularities associated 
with the onset of new subbands [32] • At these points, the concept of i?/L-movers breaks 
down, and some of our conclusions below may change. 

Next we express the Rashba hamiltonian (|24jl in terms of i?/L-movers. The first 
term results in 

H l = — £M") / 7r c L( k )^c nr (k). (35) 

The presence of the factor k±(n) results from a careful treatment of the phases in 
Eq. ([34)1 . In Eq. (|3*3j) we omit an additional term mixing right- and left-movers. This 
term contains a rapidly oscillating factor e ±2lQnX and therefore is strongly suppressed 
by momentum conservation. The second term in Eq. (J24)) again contains the oscillating 
phase factor e ±i « n±9n+1 ) a! , which leads to a drastic suppression of H^ at low energies and 
long wavelengths. Of course, this argument relies in an essential way on the smallness 
of the coupling u_l, as one expands around the hamiltonian H . We conclude that away 
from van Hove singularities, the only important Rashba term is given by i/J in Eq. (|35|) . 
This term has the appearance of a static homogeneous but channel-dependent magnetic 
Geld. 



4. Oscillatory TMR effects in nanotubes 

In this section we will analyze the consequences of our findings regarding spin-orbit 
couplings in CNTs, see Eqs. ()35|) . for the observability of spin precession effects encoded 
in the Datta-Das oscillations of the TMR. Based on our expressions, it is possible to 
estimate the order of magnitude of this effect. 

For a concrete estimate, let us put E = 0.2eVc/ (kD), where D is the gate-tube 
distance, Vq the gate voltage, and k denotes the dielectric constant of the substrate. 
For a given channel n, the Rashba-induced energy splitting is then easily estimated as 

AE n 0.6dv In + n | a 2 (3d A \l 

—T7- = (71 + 72) 5 -j n , 

eV G /j, R AkD 

where A c = H/mc = 3.86 x 10~ 13 m is the Compton length. Plugging in the definition 

of a, j3, we get 

AK = 0-6(7i+ 72) (d/ao)5J Xg_ |n + np|_ 
eVc 25Qttk Doq kpR 

Bands with small n are only weakly split, and hence do not contribute to oscillatory 
TMR behavior. This argument suggests that Datta-Das oscillations in principle could 
survive in a CNT, even when there are many channels. The major contribution will 
come just from the few bands with the largest n. 

To estimate the accumulated phase difference due to the different precession length 
of the two split eigenstates, let us put (n + noj/ikpR) — > 1, which represents the 
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dominant contribution, and set k = 1. Then Eq. (J36j) gives as order-of-magnitude 
estimate 

AE/(eV G ) « 2 x 1(T 6 a /£>. (37) 

Even when assuming a very close-by gate, this gives only a tiny splitting, in retrospect 
justifying perturbation theory. This splitting now translates into a momentum splitting 
Ak n = AE/v n , and hence into a precession phase mismatch along the CNT of length L 
[T3"j . For the nth band, this phase difference is 

A0„ = Ak n L « 2 x 10" 6 ^-^^. (38) 

D hv n /a 

This phase difference should be of order 2ir to allow for the observation of Datta-Das 
oscillatory TMR effects |13j . 

Away from a van Hove singularity, Eq. ()38j) predicts that oscillations appear on a 
gate voltage scale of the order of 10 6 to 10 7 V for L w D, which would make Datta- 
Das oscillations unobservable. This argument also shows that this interpretation can be 
ruled out for the parameters relevant for the Basel experiment JU]- From Eq. ([38)1 . we 
can then suggest several ways to improve the situation. First, one should use very long 
CNTs, while at the same time keeping the gate very close, and second, an enhancement 
can be expected close to van Hove singularities. Of course, very close to a van Hove 
singularity, some of our arguments above break down, but the general tendency can 
nevertheless be read off from Eq. (|38|). Furthermore, electron-electron interactions can 
also enhance spin-orbit effects [22| 125]. 

To conclude, we have presented a detailed microscopic derivation of Rashba spin- 
orbit coupling in carbon nanotubes. It turns out that the Rashba SO coupling is 
small, and therefore the prospects for observing spin-precession effects like Datta-Das 
oscillations in the tunneling magnetoresistance are not too favorable. However, for very 
long CNTs, close-by gates, and in the vicinity of a van Hove singularity, the requirements 
for observability of these effects could be met in practice. 
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